Channeling chaotic transport in a wave-particle experiment 
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I. INTRODUCTION 

The interaction of a beam of charged particles with 
electromagnetic waves is ubiquitous in physics, and it is 
central to many useful devices such as particle acceler- 
ators, plasma fusion experiments or free electron lasers. 
In these experimental set-ups, the waves are used to ac- 
celerate the particles or to guide them by assigning a 
specific mean trajectory. However, the dynamics of these 
systems is usually characterized by the competition of 
many degrees of freedom and thus, shows generically 
chaotic behaviors. Such behaviors represent a limit to 
the achievement of high performances in these devices. 
Consequently, there has been a growing interest in con- 
trolling chaos which here means to reduce it when and 
where it is undesirable and to increase it when it is useful. 

The sensitivity of chaotic systems to small perturba- 
tions triggered a stronginterdisciplinary effort to control 
chaos [I E 13, II H BM After the seminal work on 
optimal control by Pontryagin |?j , efficient methods were 
proposed for controlling; chaotic systems by nudging tar- 
geted trajectories E Ell El El Q • However, for many 
body experiments such methods are hopeless due to the 
high number of trajectories to deal with simultaneously. 

It was recently proposed a local control method [T^ 
which aims at building barriers in phase space and hence 
confines all the trajectories rather than following them in- 
dividually. These barriers are stable structures in phase 
space like for instance invariant tori, which are gener- 
ically broken by the perturbation. The reduction of 
chaotic behaviors is achieved by using a small apt pertur- 
bation of the system which keeps its Hamiltonian struc- 
ture. 



t Unite Mixte de Recherche (UMR 6207) du CNRS, et des univer- 
sities Aix-Marseille I, Aix-Marseille II et du Sud Toulon- Var. Labo- 
ratoire affilie a la FRUMAM (FR 2291). Laboratoire de Recherche 
Conventionne d u CEA fDSM-06-35) 
*e-mail address macor@up.univ-mrs.fr 



In this article, we consider a traveling wave tube 
(TWT) specially designed to study wave particle interac- 
tion which is used to investigate experimentally the con- 
trol method and its robustness. The dynamics in this 
experimental apparatus can be accurately represented 
using a Hamiltonian which describes the motion of a 
charged particle (with unit mass) interacting with two 
electrostatic waves : 



H(p,x,t) = — + £\ cos{k\x — uj\t + cpi) 



(i) 



where (p, x) G R x [0, L] are the momentum and position 
of the particle in a tube of length L. The amplitudes, 
wave numbers, frequencies and phases of the two waves 
are denoted respectively hi, coi and cpi for i — 1,2. 
We notice that the beam intensity is sufficiently low such 
that the wave growth rate is negligible upon the length of 
the experiment that is we are in the test-particle regime. 

Generically, the dynamics of the particles governed 
by Hamiltonian is a mixture of regular and chaotic 
behaviors, mainly depending on the amplitudes of the 
waves. The Chirikov parameter [Tjj defined as the ratio 
between the two half- width of the primary resonances by 
the distance between these resonances, i.e., 



(2) 



\uj 2 /k 2 - c^i/feil 



gives a first rough approximation of the chaoticity de- 
gree of the system. Hamiltonian has a typical behav- 
ior of integrable system for small values of this param- 
eter (s<C 1). For large enough amplitudes of the waves 
(s ~ 1), large scale chaos occurs in phase space. As a con- 
sequence, the particle can have an arbitrary velocity in 
between the two phase velocities of the waves (uo 2 / k 2 and 
oj\/k\). In this TWT, such typical chaotic behavior has 
been observed directly pj|- This chaotic diffusion of the 
particles in phase space can be reduced by using an apt 
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control term which consists here as an additional wave 
(or more generally a set of waves) of small amplitude. 
The characteristics of this additional wave are computed 
explicitly, and then the wave is injected in addition to 
the two others. The results presented in this article were 
announced in Ref. flf^ . 

The paper is organized as follows : The control method 
is briefly recalled in Sec. Ill Al and its application to the 
considered Hamiltonian is presented in Sec. Ill Bl Nu- 
merical investigations of the effect of the control term 
and its robustness are reported in Sees. Ill CI and III Dl 
In Sec. IIII Al a description of the experimental set-up 
precedes the results of the implementation of the con- 
trol term shown in Sec. IIII 51 as well as its robustness in 

Sec, mm 



II. LOCAL CONTROL METHOD APPLIED TO 
A TWO WAVE MODEL 



The Hamiltonian of an integrable system can be lo- 
cally written as a function Ho(A) of the action variables 
A = (Ai,A2,...,Ad) G i.e. it does not depend on 
the conjugate angle variables = #2, . . . , 9d) G T d = 
M d /(27rZ) d , where T d is the d-dimensional torus param- 
eterized by [0,27r[ d . The equations of motion for Ho (A) 
show that the action variables are constant, and conse- 
quently the trajectories with given actions A are con- 
fined to evolve on a <i-dimensional torus with frequency 
vector cc?o(A) = dHo/dA. The dynamics on this torus is 
periodic or quasi-periodic: 0(t) = u>o(A)t + 0(0) with 
frequency vector ujq(A). In the particular case given 
by Hamiltonian an integrable situation is given by 
S\ = £2 = so that the dynamics of the integrable 
system H = p 2 /2 is characterized by constant velocity 
(p = const). A monokinetic beam of charged particles 
remains monokinetic. 

If the system described by Hq is perturbed, i.e. we 
consider the Hamiltonian 



A. Expression of the local control term 

The aim is to provide an explicit expression for an ad- 
ditional perturbation such that a specific invariant torus 
is reconstructed in the modified system. We state here 
the main result which has been extensively described in 
Ref. 15] : We consider Hamiltonian systems written as 

H(A,0)=w-A + W(A,0), 

where u> is a non-resonant vector of R d . Without loss of 
generality, we consider a region near A = (by transla- 
tion of the actions) and, since the Hamiltonian is nearly 
integrable, the perturbation W has constant and linear 
parts in actions of order e, i.e. 

W(A, 0) = ev(0) + sw(0) • A + Q(A, 0), (3) 

where Q is of order 0(||A|| 2 ). We notice that for £ = 0, 
the Hamiltonian H has an invariant torus with frequency 
vector u? at A = for any Q not necessarily small. The 
controlled Hamiltonian we construct is 



H C (A, 0) = u ■ A + W(A, 9) + f(0). 



(4) 



The control term / we construct only depends on the 
angle variables and is given by 

f(e) = w(o,d)-w(-rd(>w(o,d),e), (5) 

where do is the derivative operator with respect to 0, 
and T is a linear operator defined as a pseudo-inverse of 
u> • <9#, i.e. acting on W = W^e lh e as 

Note that / is of order e 2 . For any function W, Hamil- 
tonian Q has an invariant torus with frequency vector 
close to uj. The equation of the torus which is restored 
by the addition of / is 

A = -Td e W(O,0), (6) 

which is of order e for W given by Eq. (0 



H(A,9) = H (A) + V(A,8), 

the integrability is generically lost and the system be- 
comes chaotic. Even if KAM theorem establishes the 
stability with respect to small perturbations of invariant 
tori with a sufficiently incommensurate frequency vector 
these tori are destroyed when the amplitude of the per- 
turbation V is large enough. The break-up of invariant 
tori leads to a loss of stability of the system until the last 
invariant torus of the integrable case is destroyed and 
then large scale diffusion occurs in phase space. In the 
case of a beam of charged particles whose dynamics is 
given by Hamiltonian QJ, for £1, £2 sufficiently large, an 
initially monokinetic beam will spread in velocity due to 
this diffusion. 



B. Computation of the control term for a two wave 
system 

We consider Hamiltonian with two waves, where 
the wavenumbers are chosen according to a dispersion 
relation k\ = K(uj\) and ki = K(uj2) plotted on Fig.^ 

In order to compute /, Hamiltonian JQ with 1.5 de- 
grees of freedom is mapped into an autonomous Hamil- 
tonian with two degrees of freedom by considering that 
t mod 2tt is an additional angle variable. We denote E 
its conjugate action. The autonomous Hamiltonian is 

p 2 

H(E,p,x,i) = E + — + s\ cos(&i£ — uj\t + ip\) 

+£ 2 cos(k 2 x - u 2 t + <p 2 ). (7) 
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This barrier of diffusion prevents a beam of particles to 
diffuse everywhere in phase space. We emphasize that 
the barrier persists for all the magnitudes of the waves 
(ei,£2). 

The control term © has four Fourier modes, 
(2fci, —2(Ji) : (2fc 2 , — 2cc; 2 ), ((fci + fc 2 ), — (^i + ^2)) and 
((fci — k 2 ),— (^1 — ^2))- If we want to restore an in- 
variant torus in between the two primary resonances ap- 
proximately located at p « w\/k\ and p « u 2 /k 2l the 
frequency cj has to be chosen between the two group ve- 
locities of the waves. If we consider a beam of particles 
with a velocity in between the velocities of the waves, 
i.e., v\ = uj\jh\ and v 2 = oo 2 /k 2l the main Fourier mode 
of the control term is 
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2(cjfci — uj\){ujk 2 — 002) 



FIG. 1: TWT dispersion relation (circles) with the helix mode 
at 30 MHz (square) and the beam mode at the same fre- 
quency but with propagation velocity choose equal to about 
2.5 x 10 6 m/s (triangle); the control wave corresponds to the 
beating of these two modes 



Then, the momentum p is shifted by uj in order to define 
a local control in the region p w 0. The Hamiltonian is 
rewritten as 

H = E + uop + £1 cos(/cix — uj it + ipi) 

^2 



x cos[(fei + k 2 )x - (ui + co 2 )t + <pi + (^2]. 

A convenient choice is a; = (ui + ^2)/2 and the control 
term is given by: 

2s 1 £ 2 



(vi - v 2 y 



h = — -^-TTo. cos[(fci + k 2 )x - (ui + o^)* + y?i + ^2]. 

(11) 

Using this approximate control term does not guarantee 
the existence of an invariant torus. However, since the 
difference between / given by Eq. and f 2 given by 
Eq. is small, it is expected that for a Chirikov pa- 
rameter s not too large, the effect of the control term is 
still effective and the barrier is restored close to 



+s 2 cos(k 2 x - UJ 2 t + (f 2 ) + 



P 



(8) 



p(x,t) 



vi + v 2 2s\ cos(/cix — uj\t + (fi) 



vi - v 2 



We rewrite Hamiltonian ((HJ into the form Q where: 

ev(x, £) = Si cos(kix — UJ\t + + £2 cos(/c2X — uo 2 t + 

^(x, t) = 0, 

Q(p,M) =//2. 

The frequency vector of the selected invariant torus is 
uj = (cj, 1). From Eq. JSJ) we have that / is given by 



2e 2 cos(/c 2 x - CJ 2 t + cp 2 ) 



^2), 



^1 - ^2 



Numerical results 



(12) 



/(x,t) = --(r^) 2 , 



which is 



f(x,t) = 7 £lfel — cos(/cix-^it + (^i) 
2 L o;/ci — cji 



^2&2 



■ COS 



(/c 2 X-^ 2 ^ + ^2)] 2 , (9) 



uk 2 — LU 2 

provided uo ^ W\jk\ and uo 7^ uo 2 /k 2 . Adding this exact 
control term to Hamiltonian QJ, the following invariant 
rotational torus is restored : 

p(Xjt) = CO COs(/uiX — UJ\t + 

— (J 1 



In this section we perform a numerical investigation 
of the effect of the exact and approximate control terms 
on the electron beam dynamics. We introduce the pa- 
rameter r given by the ratio of the two wave amplitudes 
r = e\je 2 . In order to reproduce as close as possible 
the experimental set-up described in the next section 
(see also 18J), we consider the following values of am- 
plitudes, wave numbers, frequencies and phases of the 
two electrostatic waves: (ei, fei, lji, opi) = (er, 1,0,0) and 
(£2,^2,^2,^2) = (£, fc,fc,0). Thus Hamiltonian JJ) can 
be written as 

• er cosx + ecos[k(x — t)], (13) 



H(p,x,t) = — 



e 2 k 2 



i.e. v\ = and v 2 = 1. We perform simulations with 
r = 0.082 and k = 5/3. The amplitudes of the waves are 
determined by r and e that are related to the Chirikov 
parameter by the following equation: 



ojk 2 - uj 2 



cos(k 2 x - uo 2 t + (f 2 ). (10) 



s = 2v^(V^+l). 



(14) 
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The value of e will be given by s. In the following we 
consider two values of s, that is s = 0.85 (e ~ 0.11) and 
s = 1.27 (e ~ 0.24). In this case the expression of the 
exact control term given by Eq. (J2J becomes 



f(x, t) = —2e 2 (r cos x — cos k(x — t)) z 



(15) 



while the approximate control term given by Eq. ([TT]) is 



f(x,t) = 2e 2 rcos[(k + l)x - kt]. 



(16) 



Poincare sections of Hamiltonian ([T3]) computed for two 
values of the Chirikov parameter, s = 0.85 and s = 1.27, 
are depicted in Figs. 12131 panels (a). We notice that in 
both cases no rotational invariant tori survive and there- 
fore trajectories can diffuse over the whole phase space in 
between the two primary resonances. As expected, when 
the exact control term given by Eq. JTSJ) is added to the 
original Hamiltonian, then rotational invariant tori are 
restored. This is shown by the two Poincare sections in 
Figs. El and panels (b), corresponding to the two val- 
ues s = 0.85 and s = 1.27 of the Chirikov parameter. 
We notice that s = 1.27 corresponds to a chaotic regime 
where the two resonances overlap according to Chirikov 
criterion (s > 1). Nevertheless the exact control term is 
able to reconstruct the invariant torus predicted by the 
method and to regularize a quite large region around the 
recreated invariant torus. 

In order to study the effect of a simplified control term 
on the electron beam dynamics we perform numerical 
simulations adding the control term given by Eq. ([THj) 
to Hamiltonian ([T3|). As one can see from the Poincare 
section depicted in Fig. [2J panels (c), the effect of the 
approximate control term is still present with the recre- 
ation of a set of invariant tori for s = 0.85. However, this 
regularization apparently disappears when we consider 
the fully chaotic regime with s = 1.27 (see Fig. [31 panel 
(c)). Nevertheless the approximate control term has still 
a significant effect on the reduction of chaotic diffusion. 
This fact can be observed on the probability distribution 
functions of the electron beam velocity. This diagnostic 
will also be used in the experiment in order to see the 
effect of the control terms. 

In Figs. 21 and [5J the initial velocity distribution func- 
tion of a set of 10 4 particles is compared with the final one 
obtained by integrating over a time t = 50, the dynam- 
ics governed by Hamiltonian (|13|) without control terms, 
plus the exact control term JTSJ) and plus the approxi- 
mate control term (|16|) . This investigation is performed 
for the two different values of the Chirikov parameter, 
s = 0.85 for Fig. E|and s = 1.27 for Fig. In the case 
without any control term the original kinetic coherence 
of the beam is lost which means that some electrons can 
have velocities in the whole range in between resonances, 
v G [—0.67, 0.88]. Adding the exact control term the par- 
ticles are confined in a selected region of phase space by 
the reconstructed invariant tori, and the beam recovers 
a large part of its initial kinetic coherence. For s = 0.85, 
velocities of the electrons are now between 0.33 and 0.79. 
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FIG. 2: Poincare sections of Hamiltonian (|13j) for s = 0.85 
without control term (panel (a)), plus the exact control 
term (|15|) (panel (b)) and plus the approximate control term 
O (panel (c)). 
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FIG. 3: Poincare sections of Hamiltonian (|13j) for s = 
1.27 without control term (panel (a)), plus the exact con- 
trol term (|15|) (panel (b)) and plus the approximate control 
term (|16|) (panel (c)). 



FIG. 4: (a): initial beam velocity distribution function for s — 
0.85. (b): final beam velocity distribution function without 
control term @ (light gray line), with the exact control term 
(black line) and with the approximate control term (fTS)) (dark 
gray line). 

The exact control term is also efficient in the fully chaotic 
regime (s = 1.27). Concerning the approximate control 
term it is very efficient for s = 0.85 while its efficiency is 
smaller in the strongly chaotic regime. However, it has 
still some regularizing effect, inducing the reconstruction 
of stable islands in phase space which can catch and thus 
confine a portion of the initial beam particles. 



D. Robustness of the method 

The robustness of the control method for the case 
s = 1.27 is studied with respect to an error on the phase 
or on the amplitude of the computed control term. In ex- 
periment, given the frequency uo\ +002, the wave number 
ki + &2 of the control term does not satisfy in general the 
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FIG. 5: (a): initial beam velocity distribution function for s — 
1.27. (b): final beam velocity distribution function without 
control term @ (light gray line), with the exact control term 
(black line) and with the approximate control term (fTfiJ) (dark 
gray line). 



dispersion relation k = K(uo) since the dispersion relation 
is not linear. In our case it means that the experimen- 
tally implemented control term is not the exact one. For 
this reason we investigate the robustness of the control 
term given by Eq. (|Ttj|) with a phase error (/?, that is 



f(x, t) = 2e 2 r cos[(fc + l)x -kt + ip], 



(17) 



and with an error on its amplitude ruled by a factor 5, 
that is 



f(x, t) = 2e 2 r5 cos[(k + l)x - kt\. 



(18) 



The values given by Eq. (|Tfj|) are (p = and 5=1. In 
order to quantify the robustness of the approximate con- 
trol term given by Eq. ((17)) or Eq. (|18|). we introduce 
the kinetic coherence indicator defined as the ratio of the 



FIG. 6: Kinetic coherence versus the phase introduced in the 
approximate control term given by Eq. (fT7)> . The red dash- 
dotted line indicates the kinetic coherence value for the non- 
controlled case and the dotted line the kinetic coherence value 
for the approximate control term. The overlap parameter is 
s=1.27. 



variance of the initial beam over the variance of the dis- 
tribution function after a given integration time. The 
number of particles, the integration time and the initial 
conditions are equal to the ones used in the previous sec- 
tion. 

In Fig. El we show the kinetic coherence as a func- 
tion of the phase of the approximate control term for 
the strongly chaotic regime s = 1.27. We notice that cp 
or — (p will give the same velocity distribution function for 
symmetry reason. Therefore we only consider the range 
(p G [0,7r]. The efficiency of the approximated control 
term is very sensitive with respect to the phase. In fact 
an error of 5° — 6° causes a decrease of the kinetic co- 
herence of about 50% and with an error greater than 30° 
the kinetic coherence drops in the range of values of the 
non-controlled case. 

Concerning the robustness with respect to an error on 
the amplitude of the approximate control term, we plot 
on Fig. [7| the behavior of the kinetic coherence as a func- 
tion of the 5-factor which multiplies the amplitude of 
the approximate control term. We notice that around 
the reference value of 5 = 1 (no error) there is a region 
(5 G [1, 1.3]) where the approximate control term is very 
efficient in confining the beam of test particles with a 
kinetic coherence in between [0.096,0.12]. On the other 
hand reducing the amplitude of the control term, i.e. its 
energy, there is a region where one has still a confining 
effect on the beam particle. For example for S = 0.6 
the kinetic coherence is larger by 50% of the value of the 
non-controlled case. 
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FIG. 7: Kinetic coherence versus 5 for the approximate con- 
trol term given by Eq. (|18f> . The red dash-dotted line indi- 
cates the kinetic coherence value for the non-controlled case 
and the dotted line the kinetic coherence value for the approx- 
imate control term given by Eq. (fT6)) .The overlap parameter 
is s = 1.27. 



III. EXPERIMENTAL TESTS 

A. Experimental set-up 

The experimental implementation of the control term 
is performed in a long traveling wave tube (TWT) USUI 
extensively used to mimic beam plasma interaction (22, 
12^ and recently to observe resonance overlap responsible 
for Hamiltonian chaos P3 The TWT sketched in Fig.0 
is made up of three main elements: an electron gun, a 
slow wave structure (SWS) formed by a helix with axi- 
ally movable antennas, and an electron velocity analyzer. 
The electron gun creates a beam which propagates along 
the axis of the SWS and is confined by a strong axial mag- 
netic field with a typical amplitude of 0.05 T which does 
not affect the axial motion of the electrons. The central 
part of the gun consists of the grid-cathode subassembly 
of a ceramic microwave triode and the anode is replaced 
by a Cu plate with an on-axis hole whose aperture de- 
fines the beam diameter equal to 1 mm. Beam currents, 
lb < 1mA, and maximal cathode voltages, \V C \ < 200 V, 
can be set independently; an example of typical velocity 
distribution functions is given in Figs. [3| and El (panel 
(a)). Two correction coils provide perpendicular mag- 
netic fields to control the tilt of the electron beam with 
respect to the axis of the helix. For the data shown in 
this article is chosen weak enough to ensure that the 
beam induces no wave growth and the beam electrons 
can be considered as test electrons. The SWS is long 
enough to allow nonlinear processes to develop. It con- 
sists in a wire helix that is rigidly held together by three 
threaded alumina rods and is enclosed by a glass vacuum 



tube. The pressure at the ion pumps on both ends of the 
device is 2 x 10 -9 Torr. The 4 meter long helix is made 
of a 0.3 mm diameter Be-Cu wire; its radius is equal to 
11.3 mm and its pitch to 0.8 mm. A resistive rf termi- 
nation at each end of the helix reduces reflections. The 
maximal voltage standing wave ratio is 1.2 due to residual 
end reflections and irregularities of the helix. The glass 
vacuum jacket is enclosed by an axially slotted 57.5 mm 
radius cylinder that defines the rf ground. Inside this 
cylinder but outside the vacuum jacket are four axially 
movable antennas which are capacitively coupled to the 
helix and can excite or detect helix modes in the fre- 
quency range from 5 to 95 MHz. Only the helix modes 
are launched, since empty waveguide modes can only 
propagate above 2 GHz. These modes have electric field 
components along the helix axis 23]. Launched elec- 
tromagnetic waves travel along the helix with the speed 
of light; their phase velocities, v^j, along the axis of the 
helix are smaller by approximately the tangent of the 
pitch angle, giving 2.8 x 10 6 m/s < v^j < 5.3 x 10 6 m/s. 
Waves on the beamless helix are slightly damped, with 
\kf\/\kf\ « 0.005 where k° = k 0r + ik 0i is the beamless 
complex wave number. The dispersion relation closely re- 
sembles that of a finite radius, finite temperature plasma, 
but, unlike a plasma, the helix does not introduce ap- 
preciable noise. Finally the cumulative changes of the 
electron beam distribution are measured with the veloc- 
ity analyzer, located at the end of the interaction region. 
This trochoidal analyzer |24| works on the principle that 
electrons undergo an E x B drift when passing through 
a region in which an electric field E is perpendicular to a 
magnetic field B. A small fraction (0.5%) of the electrons 
passes through a hole in the center of the front collector, 
and is slowed down by three retarding electrodes. Then 
the electrons having the correct drift energy determined 
by the potential difference on two parallel deflector plates 
are collected after passing through an off-axis hole at the 
back of the analyzer. The time averaged collected current 
is measured by means of a pico-ampermeter. Retard- 
ing potential and measured current are computer con- 
trolled, allowing an easy acquisition and treatment with 
an energy resolution lower than 0.5 eV. In the absence of 
any emitted wave, after propagating along the helix, the 
beam exhibits a sharp velocity distribution function with 
a velocity width mainly limited by the analyzer resolution 
as shown in Figs. 191 and ITU1 (panel (a)). For Fig. Et, the 
beam with radius 3mm is diffracted bypassing through 
the three grounded grids of a spreader [25( just after leav- 
ing the gun while for Fig. fTHlv. the beam radius is 1 mm 
and the spreader has beem removed for the sake of sim- 
plicity. 



B. Experimental implementation of the control 
term 

We apply an oscillating signal at the frequency of 30 
MHz on one antenna. It generates two waves: a helix 
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FIG. 8: Sketch of the Travelling Wave Tube: (1) helix, (2) 
electron gun, (3) trochoidal analyzer, (4) antenna, (5) glass 
vacuum tube, (6) slotted rf ground cylinder, and (7) magnetic 
coil. 




v [10° m/s] 



mode with a phase velocity equal to = 4.07 x 10 6 m/s, 
a beam mode with a phase velocity equal to the beam ve- 
locity Vb (in fact two modes with pulsation uo = kvb ± uob 
corresponding to the beam plasma mode with pulsation 
= (n^/raeo) 1 / 2 , Doppler shifted by the beam veloc- 
ity Vb, merging in a single mode since uob <^ uo in our 
conditions). Figures 191 and ITU1 fpanel (b)) show the mea- 
sured velocity distributions of the beam after interacting 
with these two modes over the length of the TWT for 
two different values of the Chirikov parameter. The case 
with s = 0.85 was previously investigated [ijj. The 
red square (resp. blue triangle) shows the phase veloc- 
ity v^ (resp.v&) of the helix (resp. beam) mode on the 
middle of the resonant domain determined as the trap- 
ping velocity width of the helix mode v^ ± 
(resp. Vb =b 2^/eCb^ /m ) where is the signal ampli- 
tude applied on the antenna and Ch® = 3542mV (resp. 
Cb<& = 286mV) is the real amplitude of the helix (resp. 
beam) mode. Both Ch and Cb are determined experi- 
mentally by the estimations of the coupling constant for 
helix Cy^resp. beam Cb) mode. As shown in Fig. El 
the helix mode coupling coefficient Ch is obtained by 
fitting a parabola through the measured upper bound 
velocity (circles) after the cold test beam with initial ve- 
locity equal to the wave phase velocity has been trapped 
by the wave at a given frequency over the total length of 
the TWT. As shown in Fig. El the beam mode coupling 
coefficient Cb is obtained by fitting a parabola through 
the measured upper bound velocity (circles) for a beam 
with a mean velocity very different from the helix mode 
phase velocity at the considered frequency. These two 
domains overlap and the break up of invariant KAM tori 
(or barriers to velocity diffusion) results in a large spread 
of the initially narrow beam of Figs. 151 and ITUlf panel (b)) 
over the chaotic region We now use an arbitrary 

waveform generator [2^ to launch the same signal at 30 
MHz and an additional control wave with frequency equal 
to 60 MHz, an amplitude and a phase given by Eq. ((11)). 
The beam velocity is also chosen in such a way that the 
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FIG. 9: Beam velocity distribution functions at the output of 
the TWT for s = 0.85 : (a) test beam (I b = 50 nA) with- 
out electrostatic wave, (b) with helix mode (trapping domain 
shown by red line with phase velocity marked by a square) 
and beam mode (trapping domain shown by blue line with 
phase velocity marked by a triangle) at 30 MHz , (c) with 
an additional controlling wave at 60 MHz and phase velocity 
given by grey circle and dotted line. 
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FIG. 10: Beam velocity distribution functions at the output 
of the TWT for s = 1.27 : (a) test beam (I b = 2 nA) with- 
out electrostatic wave, (b) with helix mode (trapping domain 
shown by red line with phase velocity marked by a square) 
and beam mode (trapping domain shown by blue line with 
phase velocity marked by a triangle) at 30 MHz, (c) with 
an additional controlling wave at 60 MHz and phase velocity 
given by grey circle and dotted line. 
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FIG. 11: Measure of helix mode coupling constant Ch cal- 
culated from the measured upper bound velocity of trapping 
domain (circles) at frequency / equal to (a) 30 MHz and (b) 
60 MHz. Black line parabola is obtained using the average 
value of Ch- Gray parabolas are given by the average error. 
The parabolic fit is only valid beyond the vertical segment 
indicating the wave amplitude beyond which beam trapping 
occurs over the length of the TWT. 



wave number of the helix mode at 60 MHz properly sat- 
isfies the dispersion relation function shown as circles in 
Fig. We neglect the influence of the beam mode at 
60 MHz since its amplitude is at least an order of mag- 
nitude smaller than the control amplitude as shown by 
comparing Figs. ITTk and E|f° r 30 MHz. As observed 
on Figs. 151 and ITUT panel (c)) where the grey circle indi- 
cates the phase velocity of the controlling wave, the beam 
recovers a large part of its initial kinetic coherence. For 
s = 0.85 (see Fig. EJ) the beam does not spread in ve- 
locity beyond the reconstructed KAM tori, in agreement 
with the numerical simulations of Fig. [21 For the more 
chaotic regime (see Fig. ITDb) with s = 1.27 the improve- 
ment of the kinetic coherence is still present as shown 
in Fig. ITUb . It can no more be associated with the re- 
construction of a local velocity barrier, as expected from 
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FIG. 12: Measure of coupling constant C&, at frequency equal 
to 30 MHz for a beam mode at 2.5 x 10 6 m/s. Black line 
parabola is obtained using the average value of C&, calculated 
from the upper bound velocity of the trapping domain be- 
tween the vertical segments where the beam become trapped 
over the length of the TWT and the "devil's staircase" is not 
yet evident |26|. Two independent measures are shown by 
circles and squares to give an error estimate. 



the numerical results in Fig. [3| (panel (c)). For this last 
overlap parameter an experimental exploration of the ro- 
bustness of the method will be shown in the next section. 



C. Robustness of the method 

In our experiment the control term is given by an ad- 
ditional wave whose frequency, amplitude and phase are 
computed as shown in Sec. |H] In order to quantify the 
robustness of the method we will compare the various ex- 
perimental situations to a reference one (Fig. ITUk). This 
reference is taken as the (initial) cold beam distribution 
function. An example of the distribution we were able 
to reach with control is given in Fig. ITUb . The control 
amplitude is 140 mV in agreement with 144 mV given by 
the method up to experimental errors. The phase is cho- 
sen experimentally and arbitrarily labelled 0°. The beam 
velocity is chosen equal to 2.498 x 10 6 m/s in agreement 
with 2.51 x 10 6 m/s as estimated from the dispersion re- 
lation shown in Fig. ^ 

We investigate the robustness of the control method 
with respect to variation of phase and amplitude in the 
approximate control term given by Eq. JTTJ. 

We use the kinetic coherence indicator to quantify the 
effect of the control, denned as the ratio of variance of 
the cold beam distribution function over the variance of 
the distribution function. Other indicators (integral and 
uniform distances) were used and gave similar results. 



Figure EH shows the velocity distribution functions for 
two values of the phase (—5° and 22.5°) keeping the other 
parameters constant. It shows that for a phase equal to 
—5° close to the reference value the two velocity distribu- 
tion functions are very similar, and more peaked at —5° 
than at 0°. For 22.5°, the control wave has the opposite 
effect, increasing chaos. In Fig. [Q] we show the kinetic 
coherence as a function of the phase of the approximate 
control term. It shows a narrow region around the refer- 
ence value where the control wave is the most efficient. 

In Fig. El we represent the kinetic coherence as a func- 
tion of the amplitude of the control wave. When chang- 
ing the control wave amplitude a resonance condition in a 
narrow region around the optimized amplitude is still ob- 
served. For amplitudes smaller than the reference (com- 
puted) value the effect of the control decays fast and the 
electron velocities are more widely spread than in the 
non-controlled case. Besides, for larger values, the beam 
velocity spread increases but the control term energy be- 
comes comparable to the beam mode energy changing 
radically the initial system. We have observed, due to 
beam current conservation, a lower peak at initial beam 
velocity implies that electron velocities are more widely 
spread. An enlargement of distribution around the main 
peak is shown in Figs.El a ?k an d confirms that 140 mV 
appears to be the optimum. 

Finally we check the sensitivity of the control mode 
with respect to the initial beam velocity. This corre- 
sponds introducing an error both on the wave number 
and on the amplitude of the control mode. The overlap 
parameter s depends on the phase velocity difference be- 
tween the helix and beam modes (see Eq.(EJ); for such 
a reason we also measure the non-controlled velocity dis- 
tribution function for each initial beam velocity. Fig- 
ure clearly exhibits the resonant condition expected 
at the reference value 2.51 x 10 6 m/s. We also note that, 
without control, chaos is continuously increasing as ex- 
pected since when the phase velocity difference decreases 
resonance overlap (and chaos) increases. Figure IT%1 shows 
how two beams with close initial velocities with similar 
chaotic behavior have two different responses to the same 
control term. 



IV. SUMMARY AND CONCLUSION 

Even if a modification of the perturbation of a Hamil- 
tonian system generically leads to the enhancement of 
the chaotic behavior, we have applied numerically and 
experimentally a general strategy and an explicit algo- 
rithm to design a small but apt modification of the po- 
tential which drastically reduces chaos and its attendant 
diffusion by channeling chaotic transport. The exper- 
imental results show that the method is tractable and 
robust, therefore constituting a natural way to control 
the dynamics. The robustness of the method has been 
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FIG. 13: Distribution functions for the maximum kinetic co- 
herence value (resp. minimum), blue line (resp. red line), is 
compared with the non- controlled distribution function (grey 
dashed line). The overlap parameter is s = 1.27. 
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FIG. 15: Kinetic coherence defined such as in Fig. 114k for dif- 
ferent control amplitudes and a phase 0° and a beam velocity 
2.498 x 10 6 m/s. Control becomes efficient in a narrow region 
around 140 mV. Solid line shows the level of the kinetic co- 
herence for the non-controlled chaos. Dashed line is only an 
eye-guide. The overlap parameter is s = 1.27. 



with an additional cost of energy which corresponds to 
less than 1% of the initial energy of the two- wave sys- 
tem. We stress the importance of a fine tuning of the 
parameters of the theoretically computed control term 
(e.g., amplitude, phase velocity) in order to force the ex- 
periment to operate in a more regular regime. For such a 
reason an iterative process to find some optimal experi- 
mental conditions is suggested for future improvement of 
the method. Other control terms can be used to increase 
stability (by taking into account the other Fourier modes 
of / given in Eq. (0 when experimentally feasible). The 
achievement of control and all the tests on a TWT as- 
sert the possibility to practically control a wide range of 
systems at a low additional cost of energy. 



FIG. 14: Kinetic coherence defined as the variance of the cold 
beam distribution function over the variance of distribution 
function for different control term phases, an optimized am- 
plitude 140 mV and a beam velocity 2.498 x 10 6 m/s. Control 
becomes efficient in a narrow region close to zero. Solid line 
shows the level of the kinetic coherence for the non-controlled 
case. Dashed line is only an eye-guide. The overlap parameter 
is s = 1.27. 

checked for an overlap parameter equal to s = 1.27 by 
changing phase and amplitude of the control term and 
beam velocity to check resonance condition on the helix 
dispersion relation. All these measurements have shown a 
significant region around the prescribed values for which 
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FIG. 16: Enlargement of velocity distribution functions close 
to initial beam velocity for different control amplitudes: (a) 
larger, (b) lower than experimentally optimized amplitude 140 
mV. (c) velocity distribution functions for three different con- 
trol amplitudes at fixed control phase equal to 0° . The overlap 
parameter is s — 1.27. 
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FIG. 17: Kinetic coherence versus initial beam velocity with 
optimized control term amplitude (140 mV) and phase (0°). 
(a), squares (resp. triangles) show the values for velocity dis- 
tribution function obtained with (resp. without) applied con- 
trol term, (b), Ratio between kinetic coherence measured 
with and without control for different initial beam velocity. 
The overlap parameter is s — 1.27. 
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FIG. 18: Velocity distribution functions for two close values of 
initial beam velocity (a) with and (b) without applied control 
. The overlap parameter is s — 1.27. 
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